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ABSTRACT 



Deep-water wave statistics for Colombo on the west coast of Ceylon 
have been compiled using the Sverdrup-Munk-Bretschne ider wave -hindcas t 
method applied to 12-hourly weather maps of the West Indian and South 
Atlantic Oceans for the one-year period from June 1968 through May 1969. 

Results of the wave -hindcas t analysis are presented in the form of 
monthly and annual he ight -period-direct ion frequency distributions. The 
predominant waves are in the one to three foot height range, have periods 
centered about 13-14 seconds, and arrive from westerly to south-south- 
westerly directions. Wave action is most frequent in May through 
September and least in December. 

The wind waves and local swell on the west coast of Ceylon strongly 
reflect the seasonal Monsoons. The principal source of distant swell 
for this coast is the prevailing westerly wind belt of the Southern 
Hemisphere between Argentina and the longitude of Ceylon (80°E) . The 
subtropical anticyclone in the central Indian Ocean is a relatively 
quiet source region for swell compared to the Monsoon belt and the 
prevailing westerlies. 
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I . INTRODUCTION 



Ceylon is an island country situated in the South Indian Ocean. 
Figure 1 shows its geographical location. Colombo (7°N, 80°E) is her 
capital and principal harbor. Lacking major natural resources, the 
island's economy is turning more and more to dependence on the sea. 

Plans are currently afoot for the construction of several small harbors 
for Ceylon's infant fishing industry which is concentrated mainly on 
the west coast. In the design of harbors, wave statistics are needed 
for determining breakwater orientations, heights, and construction 
materials, and for anticipating sanding problems caused by the effect 
of the harbor on littoral drift. It is believed that climatological 
wave data will also have other practical military and civil appli- 
cations in Ceylon in connection with coastal patrol, shoreline 
erosion, beach preservation, and related coastal engineering activities. 

Statistical wave data for Ceylon are not available and therefore 
this study was undertaken in part to compile climatological wave data 
for deep water off Colombo. It will be seen that the statistics pre- 
sented herein are applicable to most of the west coast of Ceylon. 

In carrying out this study various problems were met and some 
innovations were introduced in the application of the wave hindcasting 
technique used and in the compilation of the wave statistics. These 
are described in the following pages. In addition, the seasonal wave 
regime on the west coast of Ceylon is examined in relation to the 
synoptic and seasonal weather character of the West Indian and South 
Atlantic Oceans. 
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II. METHOD OF OBTAINING CLIMATOLOGICAL WAVE DATA 



A. AVAIIABLE METHODS 

Climatological wave data can be defined as a statistical compilation 
of the frequency of occurrence of waves of various height, period, and 
direction by months, seasons, or averaged over a year for a given 
location at sea or along a coast. The method chosen in this study 
for compiling wave statistics is the wave -hindcas t method, in which a 
series of historical synoptic weather maps are analyzed for their wind- 
field characteristics and a conventional technique is applied to fore- 
casting the waves produced which arrive at Ceylon. 

There are several wave - forecast ing methods in existence and a 
comparative description of each has been given by Timmie [1969]. In 
a broad sense they can be classified into two main groups, the spectral 
methods and the non-spectral empirical methods. 

Among the spectral approaches, the Pierson-Neumann- James [1953] 
or PNJ method is probably the most widely used. The PNJ method 
describes the sea state and the swell characteristics in the form of 
an energy spectrum, and enables the forecaster to obtain various height 
and period parameters of the waves generated. This method, or variants 
of it, is the one most frequently used by scientists interested in the 
mechanism of energy transfer and propagation as it gives more informa- 
tion on the energy distribution in the waves. However, it has a lesser 
degree of acceptance in engineering applications where the dominant wave 
heights and periods are mainly used. 
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Figure 1: GEOGRAPHICAL LOCATION OF CEYLON 

The continental shelf margin is shown by the broken line. 



There is also the French spectr o-angular densities method by Gelci, 
Cazale, and Vassal [1958] which is in common use by European scientists 
and engineers, particularly the French and the Dutch. It is not suited 
for manual synoptic work as it is extremely time consuming, but it has 
been claimed by these authors to give very good results using numerical 
analysis . 

The most widely used of the non-spectral methods is the Bre tschne ider 
modification of the Sverdrup-Munk [1947] wave-forecasting technique and 
referred to as the Sverdrup-Munk-Bre tschneider or SMB method [Bret- 
schneider, 1952]. This method provides a forecast in terms of a signifi- 
cant wave height and a significant wave period. The significant wave 
height is the average height of the highest one-third of the waves. The 
significant wave period is the period associated with the significant 
wave height, and closely approximates the period of maximum energy 
concentrat ion . 

Evaluations of the PNJ and SMB methods have been made by Timmie 
[1969], Rattray and Burt [1956], and Wiegel [1966]. According to 
Wiegel (p. 239), "One thing is apparent from the few comparisons of 
forecasts with observed waves: no major procedure is appreciably 

better than the others in use today." 

B. CHOICE OF SMB METHOD 

For this study the SMB wave - forecast ing technique was adopted mainly 
for the following reasons : 

(1) The method is simpler and quicker than other methods and 
appears to give forecast results that are equally as good. 
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(2) The applications of this study in Ceylon will involve mostly 
coastal engineering problems in which data on the dominant height and 
period are ordinarily adequate. 

(3) Even though the PNJ method provides a more complete des- 
cription of the sea, its format of wave -hindcast presentation in 
several studies in which it has been applied (National Marine Con- 
sultants - California Coast [1960], Oregon-Washington Coast [1961]) 

is identical with that using the SMB method; i.e., the wave statistics 
are presented in the form of frequency of occurrence of various com- 
binations of significant height, dominant period, and wave direction. 

A more detailed or complete presentation would be extremely cumbersome 
and the time required to compile and prepare the additional data would 
be significant. 

(4) This investigator anticipates that additional wave studies 
will be made in Ceylon under his direction. Therefore, it was desirable 
to gain experience in a technique that is both easy to teach and simple 
to use. 
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Figure 2 : POSSIBLE WAVE GENERATING AREAS FOR COLOMBO 



III. SELECTION OF WEATHER MAPS 



A. REGION COVERED 

Prior to the selection of the weather maps to be used in hindcasting 
it was necessary to determine the regions of origin of ocean waves that 
would arrive at Colombo. Figure 2 shows the possible wave-generating 
areas based upon great-circle propagation of swell. It may be noted 
that swell can be generated as distant as the coast of Argentina, 

7500 nautical miles from Ceylon, and that it may transit the entire 
South Atlantic Ocean. The potential wave -generat ing region with 
respect to Colombo is bounded directionally on the northwest by the 
Indian Peninsula, and on the south by wave refraction across the con- 
tinental shelf along the west coast of Ceylon. 

A study of weather-map sources indicated that to obtain complete 
coverage of the possible wave -generating areas affecting Colombo, two 
separate sets of surface weather maps were necessary. The following 
two sets were chosen: 

1. 12 -hourly weather maps (at 0000 and 1200 GliT) of the northern 
Indian Ocean down to 40°S . These maps are in Mercator projection and 
were obtained from the Fleet Weather Central at Rota, Spain. 

2. 12 -hourly weather maps (at 0000 and 1200 GMT) of the two 
quadrants of the Southern Hemisphere from 90°W through 0° to 90°E. 

These maps are in polar stereographic projection and were obtained from 
the National Weather Center (NWC) in Washington, D. C. Both sets were 
analyzed by the originating agencies. Examples of these maps are shown 
in Figures 3 and 4, 
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Figure 3: EXAMPLE OF ROTA WEATHER MAP 
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Figure 4: EXAMPLE OF NWC WEATHER MAP 



B. PERIOD COVERED 



Past work on the preparation of wave statistics (National Marine 
Consultants, [1960 and 1961]; Roper, [I960]) indicates that a coverage 
of three full years or more would be desirable in order to minimize the 
effect of yearly variations in wave occurrence and thereby produce 
better wave statistics. In this study, however, because of practical 
time limitations and the difficulty of obtaining two sets of weather 
maps preferably covering the same interval of time, weather maps 
covering the period of one complete year from June 1968 through May 
1969 were analyzed. 

C. SETS OF WAVE ANALYSES DERIVED 

It was clear at the beginning of the study that two different sets 
of wave information would have to be derived from the maps for Colombo-- 
locally generated wind waves, and swell arriving from distant wind areas. 
It was further convenient and useful to separate the swell hindcasted 
from both sets of weather maps at 20°S because of (a) basically differ- 
ent synoptic wave -generating conditions prevailing on the equatorial 
side of the subtropical high-pressure cell in the Indian Ocean compared 
with the higher latitudes, and (b) frequent lack of data in the lower 
latitudes of the Indian Ocean on the Rota maps. 

Therefore the analysis was done *nd is presented as three separate 
sets of wave statistics: 

(1) Wind Waves (Sea) at Colombo (from the Rota maps). 

(2) Local Swell from the equatotial and tropical zones (25°N to 
20°S) in the northern Indian Ocean (from the Rota maps). 
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(3) Distant Swell from beyond 20°S in the South Indian and South 
Atlantic Oceans (from the NWC maps). 

The statistics compiled on all three sets of waves are presented in 
the form of monthly and annual height-period-direct ion frequency dis- 
tributions in the tables of Appendices B and C. All heights and periods 
presented in these tables are the significant height and significant 
period. All of the wave statistics derived apply to deep water directly 
off Colombo. 
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IV. HINDCAST PROCEDURE 



A. HINDCAST OF DISTANT SWELL 
1 . The Wind Field 

The fundamental parameter needed for wave hindcasting is the 
surface wind. The NWC maps contain relatively few weather reports in 
the oceans and it was necessary to derive the surface wind field from 
the surface pressure analysis. Surface wind reports were used wherever 
they were available. Ship reports are also uncommon on the Rota maps, 
but weather reports at land stations around the northern Indian Ocean 
are comparatively numerous. 

In delineating wind areas in the Southern Hemisphere from the 
sea -level pressure field, it was assumed that the surface wind blows 
10° to 15° to the right of the isobar direction. A wind direction 
was considered to be capable of producing swell if the wind had a 

bearing of 30 C or less from the great-circle route to Colombo in the 

o 

case of straight and parallel isobars, and 45 or less in the case of 
curved isobars. 

To facilitate the acceptance or rejection of a possible wind 
area (the fetch), a transparent overlay, similar to Figure 2, was 
prepared for use with the weather maps under analysis. The overlay 
contained great circles drawn from Colombo defining the direction 
sectors south (180°- 191-1/4°), south -southwest (191-l/4--213-3/4 0 ),and 
southwest (2 l3-3/4°-236- 1/4°) . Superimposed on the great-circle curves 
was a family of orthogonal lines of distance from Colombo spaced at a 
200 naul ica 1-mi le interval. The overlay was used to determine the 
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direction of the surface wind with respect to the great-circle route to 
Colombo, the direction sector from which the waves from the wind area 
would arrive at Colombo, and also to measure the fetch and the decay 
distance. It may be noted in Figure 2 that Distant Swell arrives at 
Colombo only in the narrow direction range of south to southwest. 

The geostrophic wind speed was estimated from the isobar 
spacing using a conventional geostrophic wind scale. The surface 
wind speed was then computed using the relationship: V g = x S x C 

where - Surface wind speed in knots 

V g = Geostrophic wind speed in knots 
S = Air-mass stability factor 
C - Isobar-curvature factor. 

The quantity S was evaluated in the following manner by interpolation 
of the set of stability factors presented by Sverdrup and Munk [1951], 
which is reproduced below: 



T -T (°F) 
w a 


S tabi lity 


s 




< - 7 


Stab le 


0.55 




-7 to 0 


Indifferent 


0.60 




1 to 4 


Indifferent 


0.65 




5 to 10 


Unstable 


0.70 




11 to 15 


Unstable 


0.75 




> 15 


Unstable 


0.85 




where T^ = Water temperature 
T^ ~ Air temperature. 

Because of the sparsity of wind observations in 


the Southern 


observations of the 


air and sea temperature 


were 


seldom available 



for use in obtaining the stability factor. The factor was estimated, 
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however, by making the following assumptions It was reasoned that, 
if the air trajectory prior to entering the fetch lay along a latitude 
line or approximately so (i e , the winds were westerly), the air-sea 
temperature difference would be negligible and the air mass would be of 
indifferent or neutral stability (S - 0.65); similarly, if the wind 
trajectory before entering the fetch was directly from the south, so 
that the air-mass moved over relatively warmer water as it flowed into 
lower, warmer latitudes, the air-mass would be highly unstable (S = 0.80). 
A trajectory directly from the southwest was assumed to produce a con- 
dition of intermediate stability (S ~ 0.75). A trajectory directly 
from the northwest, in which warm air would be expected to pass over 
cooler water, was considered to produce a stable air mass (S ~ 0.60). 

In synoptic situations in which the air flow prior to entering the 
fetch curved cyclonically or ant icyc lonica lly , the air-mass stability 
was estimated in the same manner. The stability factors associated 
with wind trajectories in the Southern Ocean that can produce swell 
which would arrive in Ceylon are shown diagramat ica lly in Figure 5. 

Figure 5 also shows the isobar-curvature factor C for the 
different wind trajectories. The values have been taken from the set 
of curvature factors presented by Sverdrup and Munk [1951], which is 
reproduced below; 



Air-mass 
Stabi lity 


Great 
Cyc Ionic 


Moderate 
or Straight 


Great 

Ant icyc Ionic 


S tab le 


0.85 


1.00 


1.05 


Indifferent 


0.90 


1 00 


1 10 


Unstab le 


0.95 


1.00 


1.15 



The isobar curvature, defined b/ Bretschne ider [1952], was estimated 
visual Ly in practice. 
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Figure 5: AIR MASS STABILITY AND ISOBAR CURVATURE FACTORS 

The solid arrow represents the mean wind direction in the fetch. 

The broken lines represent possible air-flow trajectories prior to entry into the fetch. 



2. Fetch and Duration 



On each weather map the pertinent fetches were delineated and 
measured using the overlay described above. As there could be more 
than one fetch on a given map, the individual fetches were labelled and 
each fetch was followed on the subsequent weather maps until its ulti- 
mate disappearance due to dying winds. 

A fetch was treated as stationary if the fetch front did not move 
by more than 3° latitude from one map to the next. In these cases the 
wind waves at the forward edge of the fetch were hindcasted in the usual 
manner by taking both the fetch and the duration into account. 

In some cases the fetch advanced toward Ceylon with the westerly 
wind flow of the higher latitudes, and if the rate of advance of the 
fetch front exceeded 3° latitude between 12-hour ly maps, the fetch was 
assumed to be a moving fetch. A 3° latitude advance of the fetch front 
between maps represents a mean velocity of 13 knots, which is the group 
velocity of 10-second period waves. In these cases, only the wind 
duration was assumed to control the wave dimensions, and the fetch 
length had no significance. This procedure greatly reduced the time 
to make the hindcasts and was considered to be reasonably accurate for 
the purpose of this study. 

Initially, each new wind field generating wave trains was con- 
sidered to have been in existence six hours prior to its first appearance 
on a weather map (i.e. the wind at the first map time had already been 
blowing for six hours). The wind duration as of the time of the next 
and subsequent weather maps wis derived using the equivalent-duration 
method of Bretschne ider [1952] applied to the intervening 12-hour 
interval between maps. 
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3 . Minimum Swell Height 



With practical applications in mind, it was decided that waves 
arriving at Colombo having a height less than one foot would be of 
minor importance. Accordingly, wind waves in any fetch found to have 
a height of one foot or less were discarded immediately. It was 
realized, in the case of Distant Swell, that it was possible for waves 
which were greater than one foot in the fetch to be lower than one foot 
upon arrival at Colombo. 

In order to recognize and discard these cases and thereby save 
the time of going through the swe 11-hindcast procedure, the following 
table was constructed which gives the decay distance for various com- 
binations of wave height in the fetch, H_, and minimum fetch, F . , 

r min 

that will yield a swell height at Colombo of one foot: 



Height in 
Fetch 




Minimum 


Fetch , 


F . (nautical miles) 
min 




Hp (feet) 


50 


100 


200 


300 400 600 


800 


5 


525 


850 


1600 


2000 2200 3200 


4000 


10 


1800 


3000 


5000 


6000 7000 




15 


3100 


5000 








20 


4500 


8000 








25 


6000 










30 


8000 











The minimum fetch, F . , as used in the SMB method, is defined 

mm 

as the smaller of the two fetches, namely, the actual fetch and the 
fetch equivalent to the wind duration (the fetch that would permit a 
wind with a particular speed to put into the sea the same energy and 
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hence generate same size waves as would be generated by the same wind 

speed blowing for the given duration). 

An example of the use of the table is as follows. If H is 

F 

10 feet and F . dquals 100 N M, then the swell height will.be reduced 
min 1 

to one foot after travelling over a decay distance of 3000 NM. All 
decay distances greater than 3000 NM would produce waves that are less 
than one foot high at Colombo. Therefore, only if the decay distance 
is equal to or less than 3000 NM would these waves be considered for 



swell computation. 

4 . Summary of Hindcast P rocedure 

The steps used in performing a wave hindcast are shown 
sequentially in the following table. The table is essentially a 
reproduction of the wave-hindcasting data sheet prepared for this study. 



Step . 

Wave train 

Date-time 

Isobar spacing 

Mean latitude 

Geostrophic wind , Vg 

Stability correction, S 

Curvature correction, C 

Surface wind, V 

Fetch, F 

Wind duration 

Wave direction at Colombo 

Height in fetch, 

Period in fetch, T^ 

Decay distance, D 
Minimum fetch, F 

Height at end of decay, H r 
Period at end of decay, T 
Travel time, At 
Expected time of arrival 



min 



D 



Obtained from 

Weather map 
Weather map 
Weather map 
Weather map 
Geostrophic wind scale 
Figure 5 & weather map 
Figure 5 & weather map 
V s « V g x S x C 
Overlay on weather map 
Weather map 

Overlay on weather map 
Plate D-l 
Plate D-l 

Overlay on weather map 
Plate D-l 
Plate D-2 
Plate D-2 



Plate D-3 

Weather map time At 



Purpose 

Identify fetch 

Calculate geo. wind 
Calculate geo. wind 
Calculate surface wind 
Calculate surface wind 
Calculate surface wind 
Determine, , 1 T^ 
Determine H , T 
Determine HI,, T* 

Final result 
Determine H 
Determine T 
Calculate uP , T & At 
Determine H^/R? T /T 
Calculate H? F D F 
Calculate T^ 

Final result 
Final result 
Estimate arrival time 
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B. HINDCAST OF LOCAL SWELL 



The tropical and equatorial zones of the western Indian Ocean shown 
on the Rota weather maps were analyzed for Local Swell. The potential 
swell-generating region is bounded by the south and west-northwest 
directions (Figures 1 and 2). Due to the sheltering effect of the Indian 
Subcontinent, no swell of any significance arrives at Colombo from 
directions to the north of west -northwest . For the purpose of facili- 
tating the analysis, an overlay similar to that described earlier was 
prepared for the Rota weather maps showing direction sectors and lines 
of equal distance (every 150 NM) from Colombo. 

In low latitudes the calculation of the geostrophic wind velocity 
from the isobaric spacing is unreliable. Hence, the surface wind was 
obtained only from reported winds, and, in some cases when wind 
observations were lacking, by objectively comparing the isobaric 
spacing with adjacent maps in which the pressure field could be related 
to reported winds. Except for the manner of obtaining the wind speed, 
the analysis was the same as that described for Distant Swell. 

C. HINDCAST OF SEA 

As in the case of Local Swell, only observed surface winds were used 
in hindcasting wind waves at Colombo. The wind speed on most weather 
maps was fairly low, averaging about 10 knots, and usually gave rise to 
the condition of the fully arisen sea. When the seas were not fully 
arisen, which was sometimes the case with strong winds, the wind duration 
ordinarily controlled the wave dimensions. 

The deep-water location for which the wave hindcasts were prepared 
lies off the continental shelf some 15 NM seaward of Colombo; however, 
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the wind-wave hindcasts effectively apply over the shelf as well and 
do not have to be decayed to a selected coastal site. Of course, both 
wind waves and swell must be modified for the effects of shoaling and 
refraction in order to obtain their characteristics at Colombo or at 
any other selected site in shallow water. 
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V. STATISTICAL WAVE DATA 



A. METHOD OF COMPILATION 

1 . Monthly Frequency Distributions of Sea, Local Swell, and Distant 
Swell 

Separate sets of wave data were compiled for Sea, Local Swell, 
and Distant Swell. This was considered necessary and convenient, as 
discussed earlier, due to the necessity of using two sets of weather 
maps and to the relative lack of data on the Rota maps, particularly 
with regard to hindcasting Local Swell. The statistics were compiled 
so as to show the frequency of occurrence of waves by significant 
height, significant period, and direction (H, T, and if) . Tables B-l 
to B-12 in Appendix B present these statistics as monthly distributions, 
and Table B-13 shows, in a similar format, the annual distribution of 
waves due to Sea alone. 

For the most accurate computation of wave statistics, the time 
distributions of H and T for each wave train arriving at Colombo should 
•be plotted, and notations made of the direction of arrival. From such 
graphs, the duration of occurrence of waves within different increments 
of H, T, and | could be measured. This exceedingly time-consuming 
process was circumvented by the procedures described in the following 
paragraphs . 

A given fetch, for as long as it existed, yielded one hind- 
casted set of waves (H, T, and \|f) arriving at Colombo for each 12- 
hour ly weather map. Thus, for the life of the fetch, two sets of wave 
data per day were produced. In preparing the table of Sea (wind wave) 
statistics, each set of waves hindcasted was considered to have a 
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duration of 12 hours. Accordingly, for a 31-day (744-hour) month with 
complete weather coverage, one wave set would have a frequency of 



some missing weather data, the percentage was arrived at by a slight 

modification. For example, the Sea statistics for September (Table B-9) 

were based on 44 12 -hourly weather maps, whereas there are actually 60 

12-hourly intervals in the month. This means that in September, 16 

weather maps had no data. In computing the frequencies of occurrence, 

it was assumed that the percentage frequencies for the entire month 

are the same as those for the 44-interval month. Therefore each wind- 

12 

wave set would have a frequency of occurrence of — ■ x 100 = 2.3$ 

of the month. The table in Appendix A shows the number of 12-hourly 
weather maps containing wind data on which the compilation of the wave 
statistics for Sea, Local Swell, and Distant Swell was baaed d . 



sets of waves generated each day (one per weather map) arrived at 
Colombo some days later with arrival times that were sometimes chrono- 
logically reversed from their order of generation, were occasionally 
nearly simultaneous, or followed one another by several days. Each 
swell set was tabulated by its date of arrival at Colombo. The number 
of swell sets arriving on any given day varied from none to 7, and this 
number equated to 24 hours was considered to determine the duration of 
each set on that day. 

It was assumed that when there were two or more arrivals of 
wave sets of one foot or greater on a given day (from one or more wind 
areas) there would not be a "calm 11 period on that day; if there was only 
one occurrence of waves of one foot or greater, then there would also be 




month. In the case of months with 



In the case of swell, particularly the Distant Swell, the two 
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one occurrence of waves under one foot during that day, or one "calm 11 
interval of 12-hours duration; finally, if there were no occurrences 
of waves of one foot or more, then that entire day was considered 
"calm". By "calm" is meant the occurrence of waves under one foot in 
height . 



The method used to calculate the swell frequency for a given 
month is best illustrated by an example. In February (see Table B-2) 
Distant Swell was based on 52 12-hour ly weather maps, of which 10 
intervals were found to be "calm". During the remaining 42 intervals 
there were 59 swell trains from various fetches of height one foot or 

greater. Knowing that the 59 swell sets occurred in 42 12-hour inter- 

* 42 x 12 

vals, the effective duration of each swell set was — — = 8.54 hours, 

Jy 

and the frequency of occurrence in percentage of the month was 

= 1.4$. In February there was one Distant Swell set that 
arrived at Colombo with the combination H = 5-6 feet, T = 17-18 seconds, 
and i|r = SSW; hence, the frequency of occurrence of these waves is 1.4 $ 
of the month, as shown in Table B-2 for this combination. The frequency 
figure of 1.4$ means that during February (56 12-hour ly intervals), waves 
having these specifications can be expected to occur 1.4$ of the time, or 
nearly 9.4 hours (1.4$ of 56 x 12). 

2 . Monthly Frequency Distributions of Combined Waves 

The Sea, Local Swell, and Distant Swell sets for each month 
(Tables B-l to B-12) have been combined in Tables C-l to C-12 to give 
the probable monthly occurrence for all waves in deep water off Colombo. 
As each of the three sets did not contain an equal number of days of data 
(see Appendix A), it was necessary to standardize the three wave sets 
before combining them. Thus, in cases where the wave information was 
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based on a fewer number of days than are contained in a month, the 
frequency of occurrence of all available data within that particular 
wave set were proportionately increased by a factor equal to the ratio 
of the number of 12 -hourly intervals in the month to the number of 12- 
hourly intervals of observed data. This procedure implies that the 
wave statistics obtained from a partial month of weather maps are 
representative of the entire month. 

The method used to combine the tables can again be illustrated 
by an example. The derivation of Table C-2 from the three wave-data 
sets in Table B.-2 will be reviewed for the combination H = 5-6 feet, 

T = 17-18 seconds, and \|f = SSW. When all of the wave data were com- 
piled chronologically through the month according to their arrival time 
at Colombo by the method described in the previous section, the number 
of intervals of "calm" (i.e., waves under one foot) was again found to 
be 10 (this was a coincidence because in most months the combined wave 
arrivals reduced the number of "calm" intervals). Thus the number of 
12 -hourly intervals of waves in February is 46 (i.e., 56-10). 

Pertinent hindcast information for the three wave categories is 
as follows : 



Sea: 29 12-hour ly weather maps contained data which 

yielded 3 wave sets. 

Local Swell: 22 12 -hourly weather maps contained data 

which yielded 13 wave sets. 

Distant Swell: 52 12-hour ly weather maps contained data 

which yielded 59 wave sets. 



The standardizing factors required to equate the number of wave sets to 

a month of complete weather data (56 12-hour ly intervals) are ~ = 1.93, 

2 9 

~22 ~ 2-55, = 1.08, respectively. Accordingly, the number of wave sets 
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that would be expected if there were no missing weather data are 
3 x 1.93 = 5.79 Sea sets, 13 x 2.55 = 33.15 Local Swell sets, and 
59 x 1.08 = 63.72 Distant Swell sets, giving a total of 102.66 wave 
sets in February. Therefore, each of the 102.66 wave sets had a 
duration of ^ = 5.38 hours, or occurred with a frequency of 

^56^x X 12^^ = ^*80$ °f the month. The frequencies of occurrence of 
each wave set was obtained using the standardizing factors as follows: 
Sea: 1.93 x 0.8$ = 1.5$ of the month, 

Local Swell: 2.55 x 0.8$ = 2.0 $ of the month, 

Distant Swell: 1.08 x 0.8$ = 0.9$ of the month. 

The entry in Table B-2 for the combination H = 5-6 feet, T = 17-18 
seconds, and j = SSW is 2.9$, which was obtained by combining the 
single Local Swell occurrence of this combination (frequency = 2.0$) 
with the only Distant Swell occurrence of this same combination in 
February (frequency = 0.9$). The frequency value of 2.9$, of course, 
means that during February waves with the given specifications can be 
expected to occur 2.9$ of the time during the month, or for a cumulative 
period of nearly 19.5 hours. 

3 . Annual Distributions 

Table B-13 shows the frequency of occurrence of Sea compiled 
for the entire year. Each value in the table was derived by summing 
the frequencies of occurrence shown in a given H-T-^f combination for 
all 12 months of Sea and dividing by 12. In a similar manner, the 
monthly frequency distributions of the combined waves (Tables C-l to 
C-12) have been averaged to produce Table C-13, which presents the 
probable frequency distribution of all waves that can be expected to 
occur in a given year. 
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B. WAVE STATISTICS GRAPHS 



1 . General Properties 

£ 

From the statistical wave tables in the Appendices, summary graphs 
of the monthly and annual frequency distributions of height, period, and 
direction have been plotted (Figures 6-11). Figure 6 presents a histogram 
of the annual distribution of wave heights. The most frequently occurring 
waves at Colombo are in the 1-2 foot height range. This pattern is preva- 
lent most of the year. In July, August, and November, however, the most 
frequently occurring waves are in the 2-3 foot height range (Figure 7) 
and are due to local seas at Colombo. Figure 6 also shows a cumulative 
frequency distribution, from which may be read directly the percentage 
of the year when waves greater than a specified height will be exceeded. 

As an example, it may be seen that waves in excess of 3 feet occurred 
20$ of the time. The monthly wave-height distribution graphs in Figure 
7 indicate that the highest incidence of the most frequently occurring 
waves (1-2 foot height range) is in May and June, and that the least 
wave action at Colombo is experienced in December. The largest waves 
hindcasted were in the 13-14 foot range, and occurred in April, May, 
and July. 

Figure 8 presents the annual distribution of wave periods. The 
histogram may be seen to approximate a normal frequency distribution 
and to be peaked in the 13-14 second period band, represented by swell 
waves. The 4-5 second peak at the short period end of the distribution 
is due to wind waves at Colombo. This unequal bimodal distribution of 
periods may also be seen in some of the monthly histograms of Figure 9. 

It may also be noted in Figure 9 that short-period wind waves are 
negligible in December. The longest periods hindcasted were in the 23-24 
second range and occurred in July. 
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Figure 10 shows the direction of wave arrival at Colombo to be 

✓ 

largely limited to the south-southwest and southwest. The westerly 
direction is due mainly to wind waves at Colombo. The monthly graphs 
(Figure 11) show that waves from the west occur most frequently from 
May through September. The figure also shows wind waves, dominantly 
from the west, are negligible in December. 

2 . Characteristics of Sea, Local Swell, and Distant Swell 
Examination of Tables B-l to B-12 shows that Sea at Colombo is 

characterized by 2-3 foot high waves with periods of 4-5 seconds arriving 

mostly from the west. Seas are most frequent during the months of May 

through September. 

* 

Local Swell arrives at Colombo mainly as 1-2 foot waves from 
the southwest to west-northwest with periods ranging from 8-14 seconds. 

It occurs most frequently during the months of May through August. 

Distant Swell is restricted to the three directions south, 
south-southwest, and southwest. The waves are very often in the 1-2 
foot range. The periods range from 9-24 seconds, but a high percentage 
of the energy is contained in the 14-16 second period band. The 
occurrence of Distant Swell generally increases from December to May 
and decreases thereafter. 

3 . Peak Values 

The highest wave set hindcasted for the year occurred in a May 
sea generated by a westerly wind with a speed of 25 knots. The signifi- 
cant height and period in deep water off Colombo were 13.8 feet and 10.0 
seconds, respectively. The wave height also exceeded 13 feet in April 
and July. 
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Frequency of Occurrence (per cent) 




Significant Height (feet) 

Figure 6: ANNUAL DISTRIBUTION OF WAVE HEIGHT 
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Cumulative Frequency of Occurrence (per cent) 
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Figure 7: MONTHLY DISTRIBUTION OF WAVE HEIGHT 
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Figure 8: ANNUAL DISTRIBUTION OF WAVE PERIOD 
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Significant Period (seconds) 

Figure 9: MONTHLY DISTRIBUTION OF WAVE PERIOD 
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Figure 10: ANNUAL DISTRIBUTION OF WAVE DIRECTION 
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Figure 11: MONTHLY DISTRIBUTION OF WAVE DIRECTION 



The longest period was one of 23.7 seconds in July which originated 
in a 56-knot wind field off the Weddell Sea, 6500 nautical miles from 
Colombo. Those swell arrived from the south-southwest and had a signifi- 
cant height of 5.4 feet. 

The furthest distance from Colombo at which any significant wave 
train originated was 7500 nautical miles. A train generated in this area 
by a 25-knot wind arrived from the south-southwest in deep water off 
Colombo 15-1/2 days later with a significant height of a little over one 
foot and a period of 13.4 seconds. 

The most frequently occurring wave set was H = 1-2 feet, T = 12-13 
seconds, and i|r = SSW. It arrived at Colombo for a cumulative duration of 
3.3$ of the year, or 12 days. 
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VI. WAVE CLIMATOLOGY OF WESTERN CEYLON 



A. REGIONS OF WAVE GENERATION 

The spatial window for waves arriving at Colombo includes the 
western Indian Ocean and extends completely across the South Atlantic 
in the zone of westerly winds (Figure 2). Wave generation throughout 
this window occurs in three main atmospheric circulation belts: 

(1) Trade Wind (Monsoon) Belt (20°N - 20°S), 

(2) Subtropical High Pressure Belt (20°S - 40°S), 

(3) Westerly Wind Belt (40°S - 65°S). 

These wind belts are shown diagramat ica lly in Figure 12, and their 
seasonal variations are illustrated in Figures 13 and 14. Together 
they determine the annual wave regime at Colombo and along the west 
coast of Ceylon. 

B. TRADE WIND BELT 

The most distinctive feature of the Trade Wind Belt over the Indian 
Ocean is its pronounced Monsoon circulation system, which is character- 
ized by a seasonal reversal in the prevailing wind directions caused by 
differential heating between the Asian land mass and the northern 
Indian Ocean. The phenomenon of the Monsoon is brought about in the 
following manner. 

The trade winds of both hemispheres meet in the equatorial low- 
pressure belt at a zone of convergence known as the Intertropica 1 
Convergence Zone (ICZ) or the Doldrums (Figure 12). In the Indian 
Ocean area this zone migrates over a wide range of latitudes in the 
year, as may be seen in Figure 15. It reaches its most southerly 
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position (near the Equator at the longitude of Ceylon) in February and 
its most northerly position (~ 30°N) over the Himalayas in August. During 
the months of December through March the ICZ lies to the south of Ceylon, 
and the Northeast Trade Winds, referred to as the Northeast Monsoon, blows 
over Ceylon from the Bay of Bengal. At this season of the year the winds 
blow dominantly offshore along the west coast of Ceylon, and away from 
Ceylon over the northwestern Indian Ocean. During the months of June 
through September the ICZ lies to the north of Ceylon. The Southeast 
Trades at this time of the year cross the Equator and move well into 
southern Asia, and in doing so are deflected by the coriolis effect so 
as to become a southwesterly wind, the Southwest Monsoon, over the north 
Indian Ocean, Ceylon, India, and Southeast Asia. The ICZ passes over 
Ceylon on its northward migration in May and southward again in October- 
November. These months will be referred to here as the transition months. 

As may be expected, the wave statistics for Sea and Local Swell at 
Ceylon, both of which originate in the equatorial and tropical Indian 
Ocean, reflect the seasonal change in wind direction associated with the 
Monsoons. In Figures 16 and 17 are shown the frequency of occurrence of 
waves of one foot or greater due to Sea and Local Swell, respectively. 

The figures clearly reveal the influence of the onshore winds of the South- 
west Monsoon during June through September, and the dominantly offshore 
winds of the Northwest Monsoon, as well as the irregular wind directions 
of the transition months, during the remainder of the year. The isolated 
peaks in the Local Swell in February and November may be attributed to the 
fact that the wave statistics compiled are based on insufficient data for 
these months (see Appendix A). The dominant wave direction of Sea during 
the Southwest Monsoon is from the west, as may be seen in Figure 11. 
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Figure 12 : GENERAL PATTERN OF ATMOSPHERIC CIRCULATION 

(From NAVAER 00-80U-24, 1952) 
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Figure 14: GENERALIZED PATTERN OF PREVAILING 

WINDS IN JULY AND AUGUST 

(From NAVAER 00-80U-24, 1952) 
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Figure 15: AVERAGE POSITIONS OF THE INTERTROPICAL 

CONVERGENCE ZONE IN THE INDIAN OCEAN AREA 

(From NAVAER 00-80U-24, 1952) 



46 





C. SUBTROPICAL HIGH PRESSURE BELT 

The Subtropical High Pressure Belt (20°S to 40°S) is the nearest 
source region for Distant Swell arriving at Colombo. This belt is 
represented in the central Indian Ocean by a quasi-permanent high- 
pressure cell centered at a mean latitude of 30°S around which the 

r 

wind circulation is anticyclonic (counter-clockwise). The air flow 
around the north side of the anticyclone is the Southeast Trade Wind, 
and on the south side becomes a part of the Westerly Wind Belt. The 
exposure of the west coast of Ceylon to the wind direction in this 
circulation is such that waves are not generated on its north or west 
sides. In addition, the winds in the central part of the high cell 
are normally weak and variable, so that this region is likewise not a 
source region for any significant amount of swell arriving at Colombo. 

The general lack of suitable swell-generating conditions in this 
weather belt in contrast to the Westerly Wind Belt of the more southerly 
latitudes, is reflected in the Distant Swell statistics for these two 
regions. Figure 18 shows the percentage of Distant Swell arriving at 
Colombo that originated in the latitude belts of 20°S - 40°S and 
40°S - 65°S . As may be seen, only about 7$ of the swell originated 
within the Subtropical High Pressure Belt. It may thus be seen that, 
with regard to waves arriving on the west coast of Ceylon, the Sub- 
tropical High Pressure Belt is a comparatively quiet zone lying between 
two zones of significant wave noise. 

D. WESTERLY WIND BELT 

The Westerly Wind Belt (40° - 65°S) is the major source region for 
Distant Swell arriving on the west coast of Ceylon, and is essentially 
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Figure 16: MONTHLY OCCURRENCE OF SEA 

Waves of one foot height or greater. 




Figure 17: MONTHLY OCCURRENCE OF LOCAL SWELL 

Waves of one foot height or greater. 
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Figure 18: PERCENTAGE COMPOSITION OF DISTANT SWELL 

Waves of one foot height or greater. 
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Figure 19: MONTHLY OCCURRENCE OF DISTANT SWELL 

ARRIVING FROM THE WESTERLY WIND BELT 
Waves of one foot height or greater. 
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an unbroken belt extending some 8500 NM from Argentina (60° W) to the 
longitude of Ceylon (80°E) . It is a region of intense cyclonic storm 
systems which are advected eastward with the mean air flow across the 
Southern Ocean, and which are the principal source of swell emanating 
from this region. 

Figure 18 shows that approximately 93$ of the Distant Swell arriving 
at Colombo is produced in this region. This high incidence of Distant 
Swell from these southerly latitudes indicates the importance of the 
Westerly Wind Belt as a source of waves arriving on all of the northern 
and eastern coasts of the Indian Ocean. From Figure 19 it may be seen 
that there is a seasonal variation at Colombo of the occurrence of 
swell greater than one foot propagated from the Westerly Wind Belt, 
which reaches a maxima during December to May. 

The hindcast results revealed that a predominance of the Distant 
Swell originated in the eastern South Atlantic to the south and west 
of Africa. Possible reasons for this are that this region has shorter 
decay distances to Colombo than the western South Atlantic, and also 
that the Indian Ocean to the south of Ceylon produces less swell in 
the direction of Colombo because the swell are generated by south 
winds on the west side of the eastward-moving cyclonic storm systems. 



50 



VII. UTILITY OF THE WAVE STATISTICS 



This study is the only detailed climatological wave investigation 
on Ceylon that is known to the author. Although prepared for Colombo, 
the wave statistics presented in the appendices may be applied to deep- 
water locations with reasonable judgement along most of the west coast 
of Ceylon. The wave data for directions of arrival from south through 
west can be applied without modification to any deep-water location along 
the coast. Wave statistics for directions from west to north are affected 
by the shelter of the Indian Peninsula. Along the coastal area to the 
north of Colombo, the frequency of occurrence of waves from the west- 
northwest must be reduced, whereas along the coast to the south of 
Colombo these frequencies should be increased. There will also be 
additional wave energy arriving from the northwest at stations along 
the southern half of the coast that does not reach Colombo, and which 
is not contained in the statistical tables. 

In order to apply these statistics to a shallow-water site, shoaling 
and refraction modifications must be made. These procedures are given 
in "Shore Protection, Planning and Design" [Technical Report No. 4, 

U , S. Army Coastal Engineering Research Center, 1966] and in other 
standard references on shallow-water wave transformation. 

In applying the wave statistics to coastal sites in the vicinity of 
Colombo, or elsewhere along the west coast of Ceylon, the user should be 
aware of the following limitations related to their preparation: 

1. Adequate weather data for use in hindcasting the Local Swell, 
in particular, were lacking in most months, except during the period May 
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through August (see Appendix A). In September through January and in 
April the wave statistics are based on less than six days of available 
weather data per month. 

2. Simultaneously arriving wave trains were not combined using 
energy addition to provide a resultant wave of higher significant height. 
However, in the case of small waves, as is the usual situation off 
Colombo, this would not appear to constitute a source of appreciable 
error . 

3. The representativeness of the particular 12-month period 
selected for the study is not known. One or more additional years of 
wave data would be desirable to give more reliable long-term wave 
statistics . 
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APPENDIX A 



NOTES ON THE USE OF THE STATISTICAL TABLES 
Frequency Values 

The values entered in the statistical wave tables indicate the 
percentage of the month (or year in the case of an annual table) that 
the particular wave combination is likely to occur. For example, the 
frequency figure of 1.4 for the wave set H = 1-2 feet, T = 9-10 seconds, 
and 31 SSW listed under Distant Swell for January (Table B-l, Appendix 
B) means that these particular waves can be expected to occur 1.4$ of 
the time in January, or nearly 10.4 hours (1.4$ of 24 x 31). 

Wave Height and Period 

The height and period ranges listed in the table are abbreviated 
to conserve space, and actually extend from the lower class mark of 
a height or period interval up to but not including the upper class 
mark. Thus, a wave height range of 9-10 feet means 9.0 - 9.9 feet. 

All heights are significant heights and all periods are significant 
periods . 

Data Coverage 

Although 365 days of weather maps were available, weather data 
were missing on some maps; accordingly, the wave statistics were 
determined only from the maps with data available. In preparing the 
statistical tables, it was assumed that the frequencies of occurrence 
of waves of each H-T-ty combination that were hindcasted for a partial 
month are representative of the full month, so that the percentage 
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frequencies of occurrence are identical. The number of 12-hourly weather 
maps each month from which the wave hindcasts were made are listed in 
the following table. 





Full Month 


Sea 


Local Swell 


Distant Swell 


January 


62 


18 


5 


57 


February 


56 


29 


22 


52 


March 


62 


37 


23 


59 


Apri 1 


60 


40 


10 


58 


May 


62 


54 


33 


61 


June 


60 


57 


29 


60 


July 


62 


59 


43 


50 


August 


62 


62 


31 


45 


September 


60 


44 


5 


45 


October 


62 


56 


1 


50 


November 


60 


30 


11 


55 


December 


62 


22 


4 


56 
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APPENDIX B 



MONTHLY FREQUENCY DISTRIBUTION OF SEA, LOCAL SWELL, DISTANT SWELL 

ANNUAL DISTRIBUTION OF SEA 
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Table B-l. Wave Statistics (percentage occurrence) - January 
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Table B-2. Wave Statistics (percentage occurrence) - February 
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Table B-4. Wave Statistics (percentage occurrence) - April 
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Table B-5. Wave Statistics (percentage occurrence) 
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Table B-7. Wave Statistics (percentage occurrence) -July 
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Table B-9. Wave Statistics (percentage occurrence) - September 
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Table B-10. Wave Statistics (percentage occurrence) - October 
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Table B-ll. Wave Statistics (percentage occurrence) - November 
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Table B-12. Wave Statistics (percentage occurrence) - December 
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Table B-13. Annual Sea Statistics 
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APPENDIX C 



MONTHLY AND ANNUAL 
FREQUENCY DISTRIBUTIONS OF ALL WAVES 
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Table C-l. Combined Wave Statistics (percentage occurrence) - January 
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Table C-2. Combined Wave Statistics (percentage occurrence) - February 
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Table C-3. Combined Wave Statistics (percentage occurrence) - March 




69 



Table C-4. Combined Wave Statistics (percentage occurrence) - April 
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Table C-5. Combined Wave Statistics (percentage occurrence) 
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Table C-7. Combined Wave Statistics (percentage occurrence) - July 
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Table C-8. Combined Wave Statistics (percentage occurrence) - August 
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Table C-9. Combined Wave Statistics (percentage occurrence) - September 
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Combined Wave Statistics (percentage occurrence) - October 
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Table C-ll. Combined Wave Statistics (percentage occurrence) - November 
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Table C-12. Combined Wave Statistics (percentage occurrence) - December 
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Annual Combined Wave Statistics 
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Table C-13, continued 
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APPENDIX D 



WAVE HINDCASTING CURVES 
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FETCH OISTANCE, F, IN NAUTICAL MILES 

Plate D-l: SMB SEA FORECASTING GRAPH 

(From Bretschneider , 1952) 
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Plate D-2 : SMB SWELL FORECASTING GRAPH 

(From Bretschneider , 1952) 




Plate D-3 : SMB SWELL TRAVEL TIME GRAPH 

(From Bretschneider , 1952) 
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OECAY DISTANCE, NAUTICAL MILES 
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